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O ' We apply a simple statistical mechanics cluster approximation for studying clustering in the Kern 

and Frenkel model of Janus fluids. The approach is motivated by recent Monte Carlo simulations 
work on the same model revealing that the vapor coexisting with the liquid phase contains clusters 
of different sizes and shapes whose equilibrium concentrations in general depend on the interaction 
range as well as on thermodynamic parameters. The approximation hinges on a separation between 
the intra- and inter-cluster contribution to thermodynamics, where only the former is explicitly 
^ ' computed by Monte Carlo simulations. Two levels of a simple liquid theory approximation are 

I exploited for the description of the latter. In the first we use the ideal-gas expressions and obtain 

a qualitative agreement with extensive Monte Carlo bulk simulations. This can be improved to a 
' semi-quantitative agreement, by using a hard-sphere description for the cluster-cluster correlations. 

PACS numbers: 64.60.-i, 64.70.-p, 64.70.Fx, 64.60.Ak 
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Recent advances in experimental techniques for chemical sysntcsis have provided a weU defined set of different 
^ [ protocols for obtaining colloidal particles with different shapes, chemical compositions and surface patterns. In 
particular it is now possible to obtain colloidal particles with pre-defined number and distribution of solvophobic and 
solvophilic regions on their surface. These are usually referred to as patchy colloids 

The simplest example within this realm is constituted by the so-called Janus particles, where the surface is partioned 
in only two parts with even distribution of the two philicities. In spite of their apparent simplicity, Janus particles 
^ ^ , have aroused increasingly interest in the last few years both for their potential technological applications and in view 
of the rather unusual displayed self-assembly properties as compared to conventional isotropic colloidal particlesi^"— . 

A detailed study of the fluid-fluid transition for Janus fluids has recently been carried out by Monte Carlo (MC) 
simulatioii^iiS using the Kern-Frenkel pair potentiaiii. Within this model, the solvophobic and solvophilic hemispheres 
are mimicked by an attractive square- well potential and a repulsive hard-sphere potential respectively, and two spheres 
attract each other only provided that their centers are within a given distance, as dictated by the width of the well, 
^ ' and the two attractive patches on each sphere are properly aligned one another, that is lie within a predefined relative 
angular range. 

The combined features of the equal amplitude of the two philicities coupled with the specificity of the chosen 
potential types, give rise to a micellization process originating in the vapor phase that severely destabilizes the 
condensation process thus providing a re-entrant vapour coexistence curve that in the temperature-density diagram 
is skewed toward higher densities as the system is cooled to lower temperatures^. A number of additional unusual 
features were also found for the vapour phasei^, including the fact that, for the chosen width of the square-well 
potential (50% of the particle size), there appeared a predominance of particular clusters formed by single-layer 
(micelles of about 10 particle) and double-layer (vesicles of about 40 particles) always exposing the hard-sphere part 
as an external global surface, thus inhibiting the formation of a liquid phase. 

It should be emphasized that MC simulations are particularly demanding for this system in such that very low 
temperatures (of the order of 0.25 or less in reduced units) are necessary to observe these phenomena, and this is 
expected to be even more demanding for decreasing range of the interactionsiS. 

In this paper, we focus our interest to the study of the vapor phase, following a different approach, hinging on 
a strategy similar to that devised in the context of associating fiuids, where several different theories with different 
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degrees of success have been envisioned i^srJ^ 

Our approach has been inspired by the work of Tani and Hendersoi>i^, extending the Bjerrum theory for association 
in electrolytic solutionsi^ where the total partition function is factorized into a intra- and inter-cluster contribution, 
so that the original task is reduced to the computation of the partition function for clusters of increasing sizes along 
with the interaction among them. 

While the original approaclJ^ was limited by the necessity of evaluating analytically even the intra-cluster partition 
function, in addition to the inter-cluster contribution, we propose to determine the former by explicit Monte Carlo 
simulations for each cluster and the latter using physically motivated fluid theories. Within MC simulation of each 
n— cluster, we are then able to determine the energy per particle as a function of temperature and thereby compute 
the excess free energy of the isolated cluster by thermodynamic integration. 

Clearly, this approach is particularly suited to study the vapor phase as once the first n-particle clusters have been 
simulated the resulting information can be inserted in the inter-cluster theory, and this is enough to determine the 
partition function of the vapor at all thermodynamic states. We can then follow and predict the dependence of cluster 
population on thermodynamic conditions and interaction parameters. This is particularly relevant in cases in which 
spontaneous cluster formation is particularly slow, due for example to the low value of the temperature at which 
clustering takes place, a common case when the interaction range is very short. 

The paper is organized as follows: in section we describe the model, in sections IIIII and IIVI we introduce the 
cluster theory, in section|V]we describe how we determined the intra-cluster partition function. Additional results are 
then presented in section IVTl and section IVTIl is for final remarks. 

II. THE KERN AND FRENKEL MODEL 

As in the work of Sciortino ct. alM^ we used the Kern and Frenkelii patchy hard sphere model to describe the 
Janus fluid. Two spherical particles attract via a short-range square- well potential only if the line segment joining the 
centers of the two spheres intercepts a patch on the surface of the first particle and one on the surface of the other. 
In the case of a single patch per particle, the pair potential reada^i 

$(1,2) = (/)(ri2)^'(ni,n2,fi2) , (2.1) 

where 

{+00 r < a 

-e a <r <Xa (2.2) 
Xa <r 

and 

nni,n,,h2) = ( n '[t ' " ^''"^ " (2-3) 

^ ^' ^' '^''^ \ otherwise ^ ' 

where is the angular semi-amplitude of the patch. Here ni(wi) and 112(^2) arc unit vectors giving the directions 
of the center of the patch in spheres 1 and 2, respectively, with uji ~ {61, ipi) and uj2 = {O2, (P2) their corresponding 
spherical angles in an arbitrary oriented coordinate frame. Similarly, r 12(^1) is the unit vector of the separation ri2 
between the centers of the two spheres and is defined by the spherical angle fl. As usual, we have denoted with a the 
hard core diameter and X ~ 1 + A/a with A the width of the well. 

One can define the fraction of surface covered by the attractive patch as 

X = (vl/(ni,n2,fi2))i/^., = sin^ (^|^ . (2.4) 

where we have introduced (. . ~ {l/An) J dLj{. . .) as the average over the solid angle uj. 

Reduced units ksT/e {ks is the Boltzmann constant) and pa'^ will be used as a measure of the temperature and 
density in numerical data. 

III. A CLUSTER THEORY FOR JANUS PARTICLES 

Following Ref. [l6l we split the partition function in an inter- and intra-cluster contribution. Let 7V„ be the number 
of clusters formed by n particles, where n = 1, . . . , ric (n^ being the number of different clusters) and p„ = Nn/V 
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their density. We then write the total partition function as 



?tot — 



intra\^" 
n — 1 



hnter{{Nn},V,T) , (3.1) 



where the prime indicates that the sum is restricted to all possible configurations satisfying the obvious constraint of 
conserving the total number of particles A^, 



^ n7V„ = TV . (3.2) 



n=l 



Here qlf*''^^ is the "internal" partition function for a n— particle cluster and Qintcr ({^n}, V, T) is the inter-cluster 
partition function. Additional controlled thermodynamic variables are the total volume V and the temperature T. 
The constraint can be dealt with by introducing a Lagrange multiplier so that wc minimize the quantity 

lie 

In Qtot = In Qtot + (In A) ^ n7V„ . (3.3) 

n=l 

In computing the partition function p.ip we assume that the sum can be replaced by its largest dominant contribution. 
With the help of Stirling approximation A^! (N/e)^ one then obtains 

In Qtot ~ 1^ - {Nn In 7V„ - 7V„)] + In Qi„ter ■ (3.4) 

n=l 

The correct cluster distribution {iV„} is then found from the variational condition 

d 



dNr 



■InQt 



_ = . (3.5) 

{Ar„=lV„} 



This allows the calculation of the resulting free energy, /3-Ftot = ~ In Qtot, in terms of the internal reduced free energy 
densities, -Ingjj^"'^, so that 



ric rtc rtc ^ 

= [p„ In p„ - p„] + Y + E Pn 111 ^ - V 111 ■ (3.6) 



n — 1 n—1 n—1 



In the above expression, /3 = \/{kBT). 



IV. SPECIFIC MODELS 



We now consider two specific cases, where the inter-cluster interaction is not accounted for (the ideal-gas case) or 
modeled as an effective hard-sphere like interaction between clusters. 

A. Ideal gas 

The simplest possibility corresponds to consider different clusters as non-interacting ideal particles so that 

/ y \ 

Qinter ({A^n}, V^T) ^ \{[ — ] = Qidcal , (4.1) 



n=l 



where is the dc Broglie thermal wavelength associated with each n— cluster. 
Using Eqns. p.Sp and (|4.ip one easily obtains 



^intra 

= ' (4-2) 
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where p„ = Nn/V. 

The actual value of the Lagrange multiplier A can then be numerically obtained upon inverting the constraint ()3.2 



„mtra 



E«^"^ = ^-7- (4-3) 

n=l " 



Substitution of Eq. (|4.1[) into the general expression of the free energy p.6p leads toi^ 

N N 



^ = E [ft^ i'Pn^'r) -Pr] + E Pn/^/r^^^ ■ (4-4) 



n— 1 n— 1 



B. Chemical equilibrium 



The above result (|4.4p can be used to compute chemical equilibrium among different clusters. Indeed, on defining 
/i„ as the chemical potential associated to the n— th cluster, we have 

We can then impose the equilibrium condition /i„ = ?7/ii to obtain 



(piA?)' 



(4.6) 



which can be used to compute the internal free energies, given the cluster distributions. An alternative procedure, 
based on the explicit computation of the internal energy per particle within each cluster, will be discussed in Section 

CYl 

C. Connection with Wertheim association theory 

An interesting comparison can be found with Wertheim first-order association theoryi^ which is frequently used 
in this context (see e.g. Reffiol and references therein). Within this theory, the bond contribution to the Hemholtz 
free energy can be computed from a chemical equlibrium equation under the condition that only a suitable subset of 
diagrams are included in the cluster expansion and each attractive site is engaged at most in a single bond, the limit 
of a single-bond per patch in the language of the present paper. 

Consider a system formed by only monomers and dimers, that is n = 1, 2. Then from Ea. (|4.2p and condition p.2p 
limited to n = 1, 2 we can obtain a quadratic equation in the Lagrange multiplier A. The only acceptable root can 
then be substituted into Ea. (|4.2p for ?i = 1 to obtain the fraction of patches that are not bonded, that is the fraction 
of monomers 

P 1 + V 1 + 8pA 

Using numerical simulations for n = 1, 2 clusters we are able to determine the energy per particle in an n— cluster as 
a function of temperature and thereby determine the excess free energy of the isolated cluster by integration where 
A = (q^'^''''*/[qi"''''']2)(Ai/A2)^. This equation is identical to the result from Wertheim's theory (see Eq.(lO) in ReflU) 
when translated in the appropriate language. Therefore, the present formulation is equivalent to Wertheim's theory 
provided that temperatures are sufficiently low (see Ref. [l^ for further details) and the condition single-bond per 
binding site is satisfied. On the other hand, the present theory allows for an arbitrary amplitude of the patch thus 
including the possibility of multiple bonding. 

Note that while in the case of only two clusters (n = 1, 2) requires the solution of a system of 2 coupled equations 
that results into a quadratic equation for A, a general case with clusters up to the total number of clusters ric clearly 
requires the solution of a system of ric coupled equation, a task that - in general - has to be carried out numerically. 
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D. Effective hard sphere inter-cluster interaction 



While simple, the ideal gas is clearly rather unphysical even at very low densities. A more physical description 
amounts to consider all n— particle clusters as identical hard spheres with diameters (j„ and packing fractions rjn = 
(7r/6)p„fTj^. A rather precise approximate solution in this case is provided by the Boubh'k, Mansoori, Carnahan, and 
Starling expressio n^*^'^^ , but for simplicity we here only consider the case (t„ = ctq for all n-i whose thermodynamics 
is well described by the simple monodisperse Carnahan-Starling formulae^. This can be motivated by the fact that 
only a minor variation is found in the linear cluster dimensions (sec Table H] and discussion further below) and by 
the observation that instantaneous size variations of an 71— particle cluster are comparable with the variation of the 
average cluser radii for n within a few tens. It is then attempting to approximate the correlations between different 
shaped populations of clusters by a single effective one-component hard sphere system to take care of the average 
inter-cluster correlations. 

Then 

Qintcr{{Nn},V,T) = Qidcaie"^^- , (4.8) 

where Qidcai is given in Eq. (14. ip and Fes is the Carnahan-Starling^^ excess free energy 



/gj^cs ivt) ^ m (4 - 37?t) 
Nt [l-^,f ' 

where Nt = is the total number of clusters and rjt — '?« is the total cluster packing fraction. 

Following the same steps as before one obtains 



(4.9) 



y7L 



A3 



-Givt) 



(4.10) 



where we have introduced the function 



G (x) = exp 



x{8-9x + 3x^ 



(4.11) 



For the free energy one obtains from Eq. (|3.6I) 

/3-Ftot 



V 



J2 [Pn In (P„A3 ) - p„] + J2 -PnPfn'''' + 



n=l 



n=l 



V 



(4.12) 



that differs from the ideal gas counterpart Eq. (|4.4p only for the last additional term. Clearly one recovers the ideal 
gas in the limit p„ — >■ as it should. In order to find the correct solution for this system of equation it is important 
to choose the one that is continuously obtained from the solution of the ideal gas case at ctq —)■ 0. 



E. Thermodynamic quantities 

It proves convenient to express our analysis in terms of reduced partition functions Z rather than of the full partition 
functions Q used in Section Hill This can be conveniently done by the definitions 

^ _ TT „intra _ a3 7intra (a i q\ 

(Winter - _[ ]_ .3^^ <ln ■ l^.ldj 

n=l 

Given the partition function Qtot we can determine the Carnahan-Starling excess free energy 



In 



yN 



(4.14) 



the internal energy per particle 
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+ iV ^ 2^ - AT ^ 2:5 + -^u^iT) , (4.15) 



n— 1 n— 1 

where u„ is the internal energy per particle of an cluster (see Section lTV]) . We can also determine the compressibility 
factor 

I3P _l d (In Qtot) _ 1 d (In Zjntor) _ 1 + m + V? ~ 

p p dV p dV (l-Vt)^ ' 



V. COMPUTATION OF THE INTRA-CLUSTER FREE ENERGY 

The simulation were carried out following the same prescription used for the bulk fluid phases^iiS. Two kind of 
moves for each chosen particle - a random translation and a random rotation - were allowed, following standard 
recipes^ and a standard Metropolis^^ algorithm was used to compute the energy per particle of the system of n 
particles. 

Typical runs were of about 5 x 10® steps, one step consisting of n particles moves. 

We studied first the case of clusters in the neighborhood of n = 10 particles which is expected to be sufficient to 
observe the micellization process due to the single layer clusteringiS. 

To this aim we started with an initial configuration of two pentagons with particles at their vertices juxtaposed one 
above the other. The two pentagons are parallel to the x — y plane, have the z axis passing from their centers, and 
are one at 2 = +(t/2 and the other at z = —a/2. The unit vectors attached to the spheres where chosen so to connect 
the origin to the center of the given sphere. We obtained the clusters with a lower number of particles by simply 
deleting particles and obtained the clusters with a higher number of particles by adding on the z axis a particle just 
above the upper pentagon and/or just below the lower one. However the results of the simulations are independent 
of the initial configuration chosen. 

In order to compare with previous studies^, we consider the A = 0.5(t case first. 

We performed the simulations of the isolated cluster and we have explicitly tested that results coincide with the 
calculation stemming for the bulk low density Janus fluid from which we extract cluster informations by taking all 
the cluster found with the same size and averaging their properties. 

During the simulation we allow all possible moves but we do not count the configurations which arc not topologically 
connected, i.e. those configurations where it is not possible to go from one sphere to all the others through a path; 
the path being allowed or not to move from one particle 1 to particle 2 depending whether <i>(12) has value — e or not. 

At high temperatures the limiting value for the energy per particle is — e(n — l)/n. At low temperature (/cbT/ct < 
0.15) the clusters tend to freeze into certain energy minima. This can be improved by "regularizing" the angular part 
of the Kern-Frenkel potential into 

^'(ni,n2,fi2) = {tanh[?(ni • fi2 - cos 6*0)] + l}{tanh[;(-n2 • fi2 - cos 6*0)] + l}/4 . (5.1) 

and gradually increase I starting from 1/2 during the simulation up to values where there is no actual difference 
between the continuous potential and the original stepwise one. This allowed us to reach the configuration with the 
real minimum energy with a certain confidence. 

In Fig[T] we depict the relative cluster population Nn/N as a function of the reduced density pa^ in the ideal- 
gas case for n < 12 and two different temperatures ksT/e — 0.25 (top panel) and ksT/e = 0.30 (bottom panel). 
Temperature values were selected to bracket the expected critical temperature ksT/e w 0.28 on transition from a 
vapor phase mostly formed by monomers (at higher temperatures) and a vapor phase with predominant clusters (at 
lower temperatures) in the chosen range of densities 

As expected, we observe a predominance of monomers and higher order clusters at low and high density respectively. 
No significant difference is apparent for the results of the two temperatures. This is most likely due to the ideal-gas 
nature of the interacting part and can be improved by using the Carnahan-Startling fluid description, as we shall see. 

Next we consider the internal energy per particle u„ ~ (U) jn within the 71— th cluster along with the gyration radii 
defined by 



n 



(5.2) 
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with Yav = Y^^=i^j being the position of the j— th particle. Results for both internal energy and gyration 
radii for such configurations are tabulated in Table HI This provides an additional insight on the morphologies of the 
obtained clusters, in particular on the relative weak n dependence of the linear size of the obtained clusters. 

The results for u„ as a function of temperature are reported in Table HIl and can be conveniently fitted by a Gaussian 
profile 



M„(T) = a„ exp [-6„T^] + c„ , (5.3) 

where the fitting parameters a„, &„, and c„ for the n = 2, 3, . . . , 12 clusters (ui by definition) can also be found 
in Table El 

From this expression we can determine the excess free energy of the cluster jc^^.mti'a _ ^^t'cx, mtra^^^ ^yy thermody- 
namic integration 



/:«(/?) = / dxu„(l/.T) (5.4) 
Jo 

So that fn^''"" = /^X'™tia, _|_ yki, intra ^^^-^ ^^le ideal free energy contribution being 



^id.intra ^ 3 ^ (j^^ _ ^ (5^5) 



where Vq = ttCTq/G is the volume of one n— cluster and A,i = sj 2ti fi'h?- / rUn is the de Broglie thermal wavelength, and 
with the excess part given by 



^ex>tra ^ + a„ | + ^ [crf [^KW^) - ij | . (5.6) 

The intra-cluster partition function is then Z™*'''' = Wq e^"-''"'""'"' (of course Z™^''^ = wg), where wq = 
volume of one cluster. As anticipated we here choose ct„ — ctq, for all n, where uq is the only undetermined parameter 
in the theory. 



VI. ADDITIONAL RESULTS 
A. Carnahan-Starling results 



In this case the theory depends upon the average diameter of a cluster ctq. This is obtained by the requirement 
that the Carnahan-Starling results best match MC results for the bulk simulations. 

To this aim, we consider Monte Carlo results at pa^ = 0.01 on the vapor phase, for the distribution of the cluster 
sizes, with our theory. This is depicted in Fig. [2] where we compare the Carnahan-Starling approximation with the 
MC data for the distribution of cluster sizes at decreasing values of temperatures starting from fcsT/e = 0.5 which 
provides a good match with MC results for (Tq ~ 2.64cr. This value is then used in all subsequent calculations. 

It is important to remark that, in order to find the correct solution for this system of equations, it is important to 
choose the one that is continuously obtained from the solution of the ideal gas case at ctq 0. 

At lower temperatures the discrepancy with the MC data for the vapor increases. This was to be expected in view 
of the fact that the two-layer vesiscles (n— clusters with n = 40) contribution to the vapor phase, and not included in 
the present computation, becomes increasingly importantiS. The agreement could be clearly improved by allowing a 
temperature dependence of the effective cluster diameter uq, but we have chosen to keep (Tq fixed to mantain a clear 
control of the approximations involved in our approach. 

In Fig. [3] (top panel ) we show the resulting cluster distribution for the Nn/N as a function of density for a 
temperature {ksT/e — 0.27 ) at the onset of the expected critical micelle concentration^^. Unlike previous case with 
ideal gas, there is now a clear predominance of the n w 10 clusters in the whole concentration range. Additional insights 
can be obtained by plotting the monomer density pia^ versus the total concentration pa^ for decreasing temperatures, 
as reported in the bottom panels of the same Figure, where the result of the present approach is contrasted with 
bulk numerical simulations of the same quantitj*iS. This clearly shows the onset of a critical concentration where 
clusterization becomes the predominant mechanism at each temperature (this can be obtained by extrapolating the 
flat part of the curves to the vertical exis). 

In order to assess the range of reliability of our results, we have also attempted to include in the theory all clusters of 
size up to 20 particles. Fig. Ushows how the theory compares with the MC results at ksT /e — 0.4 for the distribution 
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of the cluster sizes. Note that the vertical axis spans about 8 order of magnitudes. Here we used a slightly different 
value ctq ~ 2.92(7 for the cluster diameters. Our theory nicely follows the MC data for the vapor phase up to n < 12. 
For larger clusters discrepancies begin to show up most likely due to the fact isolated clusters tend to frequently 
disaggregate during the simulation thus providing a very low acceptance ratio. As anticipated, for this larger cluster 
sizes, a full simulation of the bulk vapor phase begins to be competitive with the present methodology, and this is the 
main reason why, in the remaining of the paper, we only consider a mixture of rt— clusters with n < 12. 

As remarked, the present theory depends upon a free parameter (the average cluster diameter (To) that is computed 
by a best fit with the bulk MC simulations. 

Fig. [5]displays the sensitivity of some of the computed quantities to the choice of the average cluster diameter cto/ct. 
In particular, we have considered the compressibility factor f3P/p, the internal energy per particle u = U/N and the 
reduced free energy per particle ln((5tot)/A^- In all cases, there is a non-negligeable dependence on the cro/f value 
indicating the importance of selecting the correct effective cluster diameter. This could be improved by considering 
a distribution of cluster diameters. 

Notice that as cto increases the packing fraction of the clusters rjt quickly exceeds unit, thus limiting the possible 
range of acceptance for the cluster diameter. Similarly, in Fig. |6]we report the compressibility factor and the excess 
internal energy per particle. The excess internal energy is compared with the MC data for the vapor phased. 

B. Prediction for a different range of the square well 

So far, we have considered the case where the range of interaction (the width of the square well) A was 50% 
of the particle size a. This is the value which has been exploited in details in past MC studies of the bulk Janus 
fluid^iiS As this range decreases, typical relevant temperatures decrease and simulations become increasingly more 
demanding from the computational point of view to equilibrate. It is then no surprising that no results have been yet 
reported in the literature for these ranges. On the other hand, these are the ranges most frequently encountered in 
the experiments^^, and this is where the usufulness of our method can be assessed. 

We have then repeated the calculations for A/cr = 0.25, that is half of previous value. 

Fig. [7] reports the cluster distributions for the ideal and the Carnahan-Starling fluids (lower temperatures), and 
are the counterpart of Figs. [l]and[3] Concentrations of the n— clusters are now shifted towards higher densities with 
respect to the case with the twice as wide range, as expected. Also now the role of the 10— cluster and the 11— cluster 
is inverted respect to before. This means that lower attractive range provides, on average, smaller stable clusters, a 
results that can be understood on intuitive basis. 

We also found that thermodynamic quantities considered above are only marginally affected by the reduction of 
the width well in the considered range of densities and temperature. 

We have also considered the case of A = 0.15cr. From Fig. [5] it is apparent that the concentrations of the rt— clusters 
are once again shifted towards higher densities respect to the case with A = 0.25(t. Also now the 7— and 8— clusters 
seems to be the ones favoured at T = 0.27 in a range of densities in a neighborhood of pa^ = 0.1. This confirms the 
trend found in the case A = 0.25(7. 



VII. CONCLUSIONS 

In this paper, we have constructed a cluster theory for the vapor of Janus fluid. This is an approach that is 
complementary to previous studies based on highly demanding MC simulations^iiS^ with the aim of providing a 
detailed description of the vapor phase in view of its remarkable unusual micellization properties. 

The main idea behind the present approach is to consider the vapor phase as formed by clusters, containing an 
increasing number of particles, that are weakly interacting among each other so that simple fluid models - such as 
ideal gas or hard spheres - can be used to mimick their physical properties. The internal degrees of freedom of each 
clusters are instead obtained through a direct MC simulation of a single isolated cluster, a much simpler task as 
compared to the bulk simulation and a procedure akin to those used in the framework of simple fluidsi^ is then used 
to combine the two calculations and obtain the full description of the system. 

It is worth noticing that, in the ideal-g similar procedure has also been already implemented in micellization 

theories by several groups^ii^, and the results we obtain in the present context are quite consistent with those. 

There are two basic reasons why we expect this approach to be valuable. First because previous full bulk simulations 
showed micelles to be only weakly interacting in the vapor density range and hence a simple description for the inter- 
cluster part is expected to be sufficient. Second, because it has been observed that the vapor properties arc mostly 
dominated by particular cluster sizes corresponding to n « 10 and n sa 40 particles, so only a limited number of 
cluster sizes is necessary to obtain a complete description. 
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In the present work, we have considered clusters up to 12 particle and compared the idcal-gas description with the 
description of a gas of hard-sphere spheres, mimicking the original clusters and with an effective cluster diameter (Tq, 
using the Carnahan-Starling approximate description. The value of ctq has been obtained by a matching of the results 
for the internal energy with full bulk MC simulations. A good agreement was found at fc^T/e = 0.5 and at densities 
pu^ = 0.01 when gq ~ 2.64cr. Results from Carnahan-Starling theory is found to be far superior as compared to the 
ideal-gas description, thus emphasizing the importance of inter cluster correlations in the vapor phase. 

We also considered higher sizes clusters (of up to 20 particles) but the agreement with the simulations for the 
larger sizes becomes less satisfactory. The theory becomes less and less accurate as oscillations in the behavior of the 
concentrations of the big clusters with size appear. This may be due to the difficulty in an accurate determination 
of the internal energy of isolated big clusters. On this respect in order to be able to observe the vesicles (clusters of 
around 40 particles^) phenomenology we certainly need to include additional insights to avoid the task of the solution 
of a system of about 40 coupled equations. An additional difficulty consists in the fact that in this case the single 
diameter effective approximation used for all clusters up to 12 in the present study, will no longer be realistic, not 
even at the simplest possible level of description. Both these problems could be tackled by focussing only on clusters 
bracketing the interesting ones {n « 10 and n « 40 in the present case). 

We showed that in accord with the simulation results of Ref. |^, at temperatures around fc^T/e = 0.27 there is 
a gap of densities where the number of clusters of 11 particles (micelles) surpasses the number of any other cluster. 
This gap shrinks as we increase the temperature. 

The determined approximation to the partition function of the vapor phase of the Janus fluid can then be used to 
compute various thermodynamical quantities. 

We found reasonable quantitative agreement between the Monte Carlo data of Ref. and our theory for the excess 
internal energy of the vapor phase of the Janus fluid. We additionally computed the compressibility factor for which 
no simulation data are yet available. 

Having validated the model against numerical predictions for A = O.Scr we pursued the analysis for lower widths 
of the well, values that are closer to the experimental range of interactions^^. In view of the overall decrease in 
the attractions, characteristic critical temperatures also decrease, thus making numerical simulations increasingly 
demanding from the computational point of view. 

For the case A = 0.25cr we produced new predictions for the concentrations, the compressibility factor, and the 
internal energy per particle as a function of density. In particular we saw that as the range of the attraction diminishes 
the Janus fluid prefer to form cluster of a lower number of particles. 

Consistent results are also found for the case A = 0.15(t, a value which rather close to those used in experiments. 

An attempt to push the cluster theory to bigger cluster sizes showed that the theory becomes less and less accurate 
as oscillations in the behavior of the concentrations of the big clusters with size appear. This may be due to the 
difficulty in an accurate determination of the internal energy of isolated big clusters. On this respect in order to be 
able to observe the vesicles (clusters of around 40 particles^) phenomenology we certainly need to include additional 
insights to avoid the task of the solution of a system of about 40 coupled equations. An additional difficulty consists 
in the fact that in this case the single diameter effective approximation used for all clusters up to 12 in the present 
study, will no longer be realistic, not even at the simplest possible level of description. Both these problems could be 
tackled by focussing only on clusters bracketing the interesting ones (n w 10 and n w 40 in the present case). 

Two additional perspectives will be the subject of a future study. First the dependence on coverage x could also 
be tackled using the present approach, and this would provide invaluable information on the micellization mechanism 
for small coverage, a task that is still our of reach to direct numerical simulations. Secondly, it would be extremely 
interesting to address the issue of the reentrant phase diagram and the (possible) existence of an additional liquid- 
liquid critical point. This has been recently attempted in a very recent preprint2S,using a monomer-cluster equilibrium 
theory in the same spirit of that presented here. 
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FIG. 1: Values of N„/N, n — 1, 2, 3, . . . , 12 as a function of the density for A = 0.5a and ksT/e — 0.25 (top panels) and 
ksT/e = 0.30 (bottom panels). In both cases curves for n = 1, . . . , 6 are on the left panels and those associated with n = 7, . . . , 12 
are on the right panels. All plots have been reported on the same scale for a better comparison. Clusters associated with values 
n = 4, 5 have curves lying below the lower limit of 0.001 concentration in the case fcsT/e = 0.25. 
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FIG. 2: Comparison between the MC data (points) and our calculations using the effective one component hard sphere inter- 
cluster partition function within Carnahan-Starling approximation for ao = 2.64(j (lines), for Nn/N, 71 = 1, 2, 3, . . . , 12 as a 
function of the clusters size n at pa^ = 0.01, A = 0.5a, and various temperatures. Also shown is the ideal gas approximation 
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FIG. 3: Values of N„/N, n = 1, 2, 3, . . . , 12 as a function of the density for ksT/e — 0.27 (top panels). Clusters with n = 1, . . . , 6 
are on the left, those with n = 7, . . . , 12 on the right. The bottom panels depict the monomer concentration pia^ as a function 
of the total density pa^ for decreasing temperatures. The result of the present approach (left) is contrasted with MC simulations 
(right). All results refer to the A = 0.5a case with a cluster diameter ao = 2.64cr. 
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FIG. 4: Comparison between MC data and cluster theory using Carnahan-Starling approximation (for cro/o" = 2.92) for Nn/N, 
n = 1, 2, 3, .... 20 as a function of the clusters size n at ksT/e = 0.4, pa^ = 0.01, and A = 0.5cr. 
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FIG. 5: Values for the compressibility factor, the internal energy per particle, and the logarithm of the total partition function 
as a function of the cluster diameter a^ja at pa'^ = 0.01, ksT/e — 0.5, and A — 0.5ct. 
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FIG. 6: Compressibility factor as predicted by the Carnahan-Starling (cro = 2.64(t) cluster theory (top panel). On the bottom 
panel we compare the MC data and the Carnahan-Starling cluster theory (same diameter as above) for the excess internal 
energy per particle for three different values of temperatures. In all cases A — 0.5a. 



17 




FIG. 7: Cluster distribution in the case A — 0.25a. The top panels report the ideal gas result at ksT/e — 0.25 {n — 1, . . . , 6 
on the left and n = 7, . . . , 12 on the right). This is the same as Fig. [T] The bottom panels depict with the same distribution 
of curves the results obtained with Carnahan-Starling approximation at fcsT/e = 0.27 which is the counterpart of Fig. [3] 
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n 


{U)/n 


U 


Rg 


1 











2 


-0.5 


-1 


-1/2 


3 


-1 


-3 ' 


-1/V3 


4 


-1.5 


-6 


0.83 


5 


-2.0 


-10 


0.76 


6 


-2.50 


-15 


0.75 


7 


-2.71 


-19 


0.91 


8 


-2.88 


-23 


0.93 


9 


-3.10 


-28 


0.96 


10 


-3.20 


-32 


1.00 


11 


-3.36 


-37 


1.04 


12 


-3.42 


-41 


1.08 



TABLE I: The low temperatures internal energy per particle of the clusters with up to 12 particles when A = 0.5cr. Also shown 
is the gyration radius Rg defined in Ea. (|5.2p . 
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-0.50 
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-0.337 


3.880 


-0.339 
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-0.346 


10.780 


-0.67 
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-0.778 


4.670 


-0.771 
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502 


-0.774 


7.975 


-0.75 
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-1.226 


5.162 


-1.025 
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890 


-1.034 


9.366 


-0.80 
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-1.700 


5.600 


-1.381 


7 


361 


-1.207 


9.214 


-0.83 


7 


-1.899 


5.263 


-1.423 


6 


767 


-1.480 


8.277 


-0.86 
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-2.064 


5.080 


-1.520 


4 


179 


-1.551 


8.503 


-0.88 
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-2.301 


5.478 


-1.579 


4 


367 


-1.681 


10.160 


-0.89 


10 


-2.394 


5.509 


-1.725 


4 


271 


-1.551 


9.419 


-0.90 


11 


-2.556 


5.644 


-1.846 


4 


829 


-1.696 


9.755 


-0.91 


12 


-2.598 


6.077 


-1.854 


5 


723 


-1.814 


10.567 


-0.92 



TABLE II: Fit to a Gaussian of the energy per particle as a function of the temperature (see Eq. (|5.3p ) 

. Cn values are common to the three cases. 



